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Abstract 

In this paper we use an 0(./V)-invariant scalar field of unbroken symmetry to investigate 
whether an interacting quantum field at the next-to-leading order Large N approximation 
may show signs of thermalization. We develop the closed time-path (CTP) two-particle 
irreducible (2PI) effective action in powers of 1/N, retaining up to next to leading order 
(0(1)) terms, and write down the corresponding (truncated) Schwinger-Dyson equations 
for its two point function. We show that in this approximation, the only translation 
invariant solutions to the Schwinger - Dyson equations are thermal. This provides a 
useful temperature concept without invoking a heat bath. Coupled with the familiar 
Kadanoff-Baym approach to quantum kinetic theory our result shows that at this order of 
approximation thermalization can occur, at least if initial conditions are smooth enough 
that a derivative expansion is valid. Our analytic result provides support for similar claims 
in recent literature based on numerical evidence. 

1 Introduction and Summary 

The problem of thermalization in relativistic quantum fields has drawn much attention over 
time, both in its own right in our attempt to understand the origin of macroscopic irreversible 
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behavior from microscopic theories, and in the context of nonequilibrium quantum field pro- 
cesses in the early Universe and in relativistic heavy ion collisions |], ||, f| ^ |5j, || [7]. 

The goal of this paper is to investigate whether at the next-to-leading order (NLO) Large 
N approximation || |], [K| [TTJ] an interacting quantum field may show signs of thermalization. 
We consider a O(N) invariant scalar field of unbroken symmetry, develop the closed time-path 
(CTP) 2-particle irreducible (2PI) p| 
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T|, m [16| effective action (EA) (TJ in powers of 
1/N, retaining up to the next-to-leading order (0(1)) terms, and write down the correspond- 
ing (truncated) Schwinger-Dyson equations for its two point function. We show that in this 
approximation, the only translation invariant solutions to the Schwinger - Dyson equations 
are thermal. Thus, without having it coupled to a heat bath this provides a useful tempera- 
ture concept. Together with the familiar Kadanoff-Baym approach to quantum kinetic theory 
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this result shows that at this order of approximation thermalization can occur, 
at least if initial conditions are smooth enough that a derivative expansion is valid. Our ana- 
lytic result provides support for similar claims in recent literature based on numerical evidence 
pT| , p2|j . We will not address the bigger question of whether the thermal solutions are in any 
sense an attractor with a nontrivial basin of attraction, as it is considerably more difficult. (We 
may look at the theory of glasses to see just how involved the long time behavior of a field 
theory can be f23|.) 

However, to appreciate this result, we must distinguish between a 'true' NLO approximation 
(truthful to the scheme) to the full quantum theory, and the theory which results from solving 
the NLO Schwinger-Dyson equations as if they were exact. By a 'true' NLO approximation we 
mean that, after deriving the Schwinger-Dyson equations up to NLO corrections, the solution to 
these equations is expanded as well, and terms of higher order which then result are discarded. 
In the second - call it the 'exact'- procedure, once the equations are written down, a solution 
is sought, which will involve terms at all orders in 1/N. 

To give an example, the situation is similar to the usual textbook derivation of the running 
of coupling constants in an interacting field theory. A renormalization group equation, let us 
say for a A0 4 scalar field theory, is derived within perturbation theory as 



dX 



eX + aX' 



(1) 



where /i is a renormalization parameter. We then solve this equation as if it were exact and get 



A(/i) 




1 



- 1 



(2) 



Call this the 'exact' way. Of course, a Taylor expansion of A (/x) should include all powers of 
A (/^o) , but only up to two powers were considered in Eq. ([]]). A true second order approxima- 
tion would be restricted to 



A (/i) 




AGuo) 1 + 



aX (/^c 




(3) 
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Call this the 'true' way. We would obtain this same solution if A 2 in the second term in Eq. 
(H) is approximated by 



A 



2 



(4) 



However, we usually believe that Eq. @ is meaningful, although it oversteps the bounds within 
which Eq. ([]]) was derived, because it can be shown that it captures the leading terms in the 
perturbative expansion (the so-called leading logs), and because in any case we only use it in 
the regime in which the coupling constant is small. 

In the thermalization problem, this procedure breaks down. If we adhere strictly to the 
'true' NLO approximation, then the only thermalization mechanism left in the theory is binary 
scattering of on-shell particles. Since this process conserves particle number, the truncated 
theory allows thermalization with non vanishing chemical potential. The 'exact' theory, repre- 
sented by the full hierarchy, on the other hand, does not conserve particle number |24j], and so 
the chemical potential must vanish in true equilibrium states. We must conclude that a "true" 
NLO approximation, in the above sense, fails to describe thermalization. 

This problem actually disappears if we solve the NLO Schwinger-Dyson equations 'exactly' 
to all orders in 1/N (like going from Eq. ([!]) to (|2|)). In this procedure the density of states 
takes on a Breit-Wigner form, and we have states with all masses. In particular, now a state 
with squared momentum — p 2 > 9M 2 [we use signature (-,+,+,+) for the Minkowski metric 
and M 2 is the leading order physical mass] may decay into three on-shell particles whereby 
particle number is no longer conserved. 

In this way, we see that the theory built on the NLO Schwinger-Dyson equations is able 
to describe thermalization, including relaxation of the chemical potential. This relaxation is a 
higher than NLO effect, and we must raise the issue of whether our analysis is still a meaningful 
approximation to the full theory. The problem is that in this case there is no analog of the 
"leading log" concept that validates the running coupling constant Eq. (^]) over Eq. (|3|). On 
the contrary, the NLO approximation discards 2^4 scattering of on-shell particles, which 
also violates particle number conservation. The contribution to the relaxation rate from these 



higher order effects is comparable to the decay of off-shell excitations p0 |. 

We therefore conclude that, while the next to leading order approximation describes ther- 
malization, it overestimates the relaxation time, and may not be a realistic picture of thermal- 
ization in actual physical systems. 



1.1 The meaning of thermalization 

Let us begin with a discussion of the exact meaning of thermalization. In the strictest sense 
an isolated system depicted by quantum field theory undergoes unitary evolution and does not 
thermalize. However, one can still ask meaningful questions such as whether certain correla- 
tion functions may converge to their thermal forms in some well defined physical limit (weak 
thermalization). 
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As a matter of fact even asking this type of questions is too ambitious. The Schwinger - 
Dyson equations for the correlation functions form an infinite hierarchy. Physical limitations 
on the precision of our measurements amounts to truncating the hierarchy. When certain 
causal boundary condition such as Boltzmann's molecular chaos hypothesis is imposed, the 
truncated subsystem will show signs of irreversibility and a tendency to equilibrate. In practical 
terms we have to deal with a truncated hierarchy to make the analysis possible. Moreover, 
nontrivial (point) field theories are plagued by divergences which can only be controlled by 
regularization and renormalization within some perturbative scheme. Therefore the question of 
weak thermalization can only make sense within a chosen approximation scheme, whether this 
is a 1/N expansion, loop expansion, expansion in powers of a coupling constant, etc., to allow 
us to organize the Schwinger - Dyson equation and evaluate the relative weight of different 
processes. 

One must distinguish between two different viewpoints. If one accepts thermalization as 
an empirical fact, then there is only the question of coming up with a formulation which 
describes this process. For example, one may assume from the beginning that the relevant 
dynamics involves only the longest scales in space and time [G5[ |2B|, |2S|, introduce a Wigner 



function as a partial Fourier transform of the propagators (as in the Kadanoff - Baym approach 



T7| , HH [T3 19) an d deduce a Boltzmann equation for the dynamics of this Wigner function 



30| , Of course the resulting model describes thermalization. 
However, this is not a proof that the original field theory thermalizes, even weakly, because 
there is no clear cut description of the set of initial conditions for which the Kadanoff-Baym 
approach, which depends upon a derivative expansion of the propagators, is valid. Mrowczynsky 
has shown that, for free theories, the only correlators that satisfy a reasonable almost-invariant 



condition are those which are exactly translation invariant f32|| . In other words, the set of 



initial conditions which permits thermalization as described by the Kadanoff-Baym equations 
may be empty, but for the equilibrium solution itself. We may say that the Kadanoff - Baym 
formalism is useful for studying certain important questions, such as the determination of the 
transport coefficients |2(| , assuming one already knows on independent grounds that the system 
thermalizes. It describes certain important aspects of the thermalization process, but it does 
not address the conditions conducive to it. 

On the other hand, if one does not assume thermalization, then one needs an approximation 
scheme which avoids the imposition of an arrow of time on the system by hand. 



1.2 The Large N Approximation Beyond Leading Order 

The number N of replicas of essentially identical fields (like the N scalar fields in an O(N) 
invariant theory, or the N 2 — 1 gauge fields in a SU(N) invariant non-abelian gauge theory) 
suggests using 1/N as a natural small parameter, with a well defined physical meaning and 
that, unlike coupling constants, is not subjected to renormalization or radiative corrections. By 
ordering the perturbative expansion in powers of this small parameter, several nonperturbative 
effects (in terms of coupling constants) may be systematically investigated. 
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In the case of the O(N) invariant theory, in the presence of a nonzero background field 
(or an external gravitational or electromagnetic field interacting with the scalar field) we may 
distinguish the longitudinal quantum fluctuations in the direction of the background field, in 
field space, from the N — 1 transverse (Goldstone or pion) fluctuations perpendicular to it. To 
first order in 1/N, the longitudinal fluctuations drop out of the formalism, so we effectively 
are treating the background field as classical. Likewise, quantum fluctuations of the external 
field are overpowered by the fluctuations of the N scalar field. In this way, the 1/N framework 
provides a quantitative measure and concrete meaning to the semiclassical approximation . 

To leading order (LO), the theory reduces to iV — 1 linear fields with a time dependent 
mass, which depends on the background field and on the linear fields themselves through a 
gap equation local in time. This depiction of the dynamics agrees both with the Gaussian 
approximation for the density matrix [[34], and with the Hartree approximation |36] . 

The ability of the 1/N framework to address the nonperturbative aspects of quantum field 
dynamics has motivated a detailed study of the properties of these systems. In non-equilibrium 
situations, this formalism has been applied to the dynamics of symmetry breaking |TT] , |37| , |38| 
and self-consistent semiclassical cosmological models f|D], [II], |4"5[ |. 

The LO 1/N theory is Hamiltonian |3(J and time- reversal invariant. However, it does not 
thermalize. For example, if we set up conditions where both the background field and the 
self-consistent mass are space-time independent, then the particle numbers for each fluctuation 



mode will be conserved. The existence of these conservation laws precludes thermalization [[44 



We note that the failure of the LO approximation to describe thermalization is indicative 
of a more general breakdown of the approximation at later times, where effects of particle 
interaction dominate. Both the distribution of energy among the field modes and the phase 
relationships (or lack thereof) among them affect the way quantum fluctuations react on the 
background or external fields. Therefore, from physical considerations, one can say that a 
theory which does not describe thermalization becomes unreliable for most other purposes as 

weii m. 

This is where the next to leading order (NLO) approximation enters. It has been applied to 
quantum mechanics [47, 48 1, classical field theory [^, |50], ^l|, |52| and quantum field theory in 1 



space dimension [M] , being contrasted both to exact numerical simulations of these systems, as 
well as against other approximations purporting to go beyond LO. The NLO has been shown 
to be an accurate approximation, even at moderate values of N. 

The 2PI formalism is also suitable for this question because, provided an auxiliary field is 
cleverly introduced, the 2PI CTP effective action can be found in closed form at each order in 

i/N m m. 



1.3 This paper 

This paper is organized as follows: In Section II we present our model, calculate the 2PI CTP 
EA, and discuss several properties of the propagators which hold to all orders in 1/N. In Section 
III we discuss translation invariant solutions to the Schwinger- Dyson equations, still without 
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any explicit approximations. In Section IV we implement the large N approximation, showing 
that it is possible to write down a closed expression for the 2PI CTP EA to next to leading order. 
In Section V we show that the only translation invariant solutions to the NLO Schwinger-Dyson 
equations are thermal. In Section VI we discuss the relaxation of the chemical potential. In the 
final Section, we ask a central question at the foundations of statistical mechanics: where does 
macroscopic irreversibility arise from microscopic reversible dynamics? We point out the exact 
spot where coarse-graining was introduced which leads to the appearance of thermalization. 



2 The model 

Let us consider a 0(iV)-invariant quantum scalar field in the limit iV — > oo. The classical 
action 

S = Jd d x (-1) j^S^S + MI^b + ^ (5) 

where <fi B) M B and A# are the bare wave function, mass parameter and coupling constant, 
soon to be renormalized, and the dimension d — 4 — e. We introduce the bare wave function 



renormalization Zb by rescaling <f)% = J(N/ZB)4> a 



s = £ fd d x {—) \d^r + Mirr + ^§- b (0> a ) 2 } (6) 

Since later on we shall discuss in detail the conservation laws of the exact and approximated 
dynamics, let us observe that this theory conserves, besides energy - momentum, a number of 
Noether charges associated with the global O(N) symmetry. Concretely, if the infinitesimal 
O(N) transformation reads 4> a — * 4> a + EaT Aol ^ then the Noether charges are 

/V r 

Q A = — d^x T Q W (7) 



Z B 

To investigate thermalization in total generality, we should allow for a Lagrange multiplier for 
each of these charges. Here we assume that all of these vanish, as well as the mean value of 
the charges themselves. Also, since there is no particle current the field is its own antiparticle, 
and the chemical potential must be zero in equilibrium. Discarding a constant term, we may 
rewrite the classical action as 



S = ^-Jd d x (-1) dpFW + 




(8) 



To set up the 1/N resummation scheme, it is customary to introduce the auxiliary field x 
writing 
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s = 



N 



Z B 
whence 



d d x 



I, 



A? 



(Ml - X ) + 



S=^-Jd d x l-ld^dPF-x 



+ 



J B 2 



2A; 



Z B 2 
(9) 



(10) 



From now on, we consider x an d 4> a as fundamental fields on equal footing. We assume a 
background field decomposition <fi = <p + (p and that the background field is identically zero 
(because of the 0(N) symmetry, the symmetric point must be a solution of the equations of 
motion) so we can focus on the dynamics of the fluctuation fields (p. We also split the auxiliary 
field x field into a background x an d a fluctuation x, x = X + X- The action becomes 



•S* — Sb ac k + S^in + S qua d + (11) 

'S'feacfc is just the classical action evaluated at cp a = 0, x — X 

9 - * 

'-'back 

S& n contains terms linear on x an d can be set to zero by a choice of the background field % 

N 



J d d x Ihf-Mlx 



(12) 



Sun = — J d d x {x-M B }x 



Squad contains the quadratic terms and yields the tree - level inverse propagators 



Squad 



N 



J B 



d a x —d^dPip* 



X a a , Zb ~2 

2** + 2~Xb" 



-x 



Finally S cu b contains the bare vertex 



Scub 



-N 



2Z 



b 



d d x {x<p a <p a } 



(13) 



(14) 



(15) 



To write the 2PI CTP EA we double the degrees of freedom, incorporating a branch label 
a = 1,2. We also introduce propagators G aa,l3b for the path ordered expectation values 



G aa ^ b (x,y) = (<p aa (x) <p?» („)) 



(16) 



and F ab for 



F ab (x,y) = ( X a (x) X b (y) 



(17) 
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Because of symmetry, it is not necessary to introduce a mixed propagator, i^x a (x) ip^ b (y)^ = 0. 
The 2PI CTP EA reads 



''back 



X 



-s, 



back 



X 



+ l -Jd d xd d y {D aa ^ b (x,y)G aa ^\x,y) + ^c ab 5(x,y)F ab (x,y)} 



~[Tr InG + Tr InF] +T C 



where c n = -c 22 = 1, c 12 = c 2 i = 0, 

D. 



N 



Cabdl - c abc x c S(x,y) 



18) 



(19) 



c a bc = 1 when all entries are 1, c a b c = — 1 when all entries are 2, and c a b c = otherwise. Tq is 
the sum of all IP I vacuum bubbles with cubic vertices from S cub and propagators G aa ^ b and 
F ab (x, y) . Observe that Tq is independent of x c - 

Taking variations of the 2PI CTP EA and identifying x 1 — X 2 — Xi we fi 11 ^ the equations 
of motion 



N 
2Zl. 



SafSCabD (x, y) - 



N 



Cab$(x,y) ~ 



ih 
~2 

ih 



G x ] aa pb (x, y) + ^U aa ^ b (x, y) = 



2X B 

N 

where D(x,y) = [d 2 x - x(x)]5(x,y), 



F- 1 ab (x,y) + ^U ab (x,y)=0 



y- {x (x) 



Ml) 



N 
2Zl 



S aP G al ^(x,x) = 



n aa ,/3b (x,y) = 2- 



sr. 



Q 



5G aa >P b (x,y)' 
We shall seek a solution with the structure 



n a b (x, y) = 2 



5F ab (x,y) 



h 



G a ^\x,y) = -5 a pG ab (x,y) 
which is consistent with vanishing Noether charges. Then it is convenient to write 



(20) 
(21) 
(22) 

(23) 

(24) 



F ab (x, y) = ^H ab (x, y) ; U aa ^ b (x, y) = 5 aP P ab (x, y) ; Ii ab (x, y) = NQ ab (x, y) (25) 



The equations become 
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^-c ab D (x, y)-i\G x l (x, y) + ^-P a b (x, y) = 



A, 



-c afe 5(x,?/) - i 



ab 



06 



AT 



(x,y) + Q ab (x,y) = 



Observe that 



±{x(x)-M*}-^-G 11 (x,x) = 



Pi \ - 2 sr Q n ( \ - 2 sr Q 

^ ab [x, y)- h 6Gab ^ y) , LJ ab [x, y)- h mab ^ y) 



(26) 
(27) 
(28) 

(29) 



These are the exact equations we must solve. The succesive 1/N approximations amount to 
different constitutive relations expressing P ab and Q a b in terms of the propagators. 

2.1 The retarded propagator 

For later use, we want an equation for the retarded propagator G ret = % (G 11 — G 12 ) = % (G 21 — G 22 ). 

(30) 



Rewrite the Schwinger - Dyson equation above as 

±-DG ac (x, V) + ^J d d z PI (x, z) G bc (z, y) = tc ac 5(x, y) 
Subtracting the (11) from the (12) components in the above equation, we obtain 



1 If 

—DG ret (x, y) + -y d d z P ret (x, z) G ret (z, y) = (-1) S(x, y) 



(31) 



where P ret = P n + Pi 



12- 



2.2 Some nonperturbative identities 

We note some non-perturbative properties of the self-energy P a b, as follows. From the identity 



5 2 T 



ipi 



ih 



G 



-1 



.. - (32) 

relating the inverse propagator to the CTP IPI EA T 1PI , the inverse propagators may be read 
off the Schwinger - Dyson equations, and T 1PI may be written as 



ipi 



N r 

— J d d xd d y 



+i 



ip al - v?" 2 



V? al - v? a2 



(x)N(x,y) 



(x) 



—D(x,y) + P(x,y) 



y? al - ^ a2 



(y)} 



tp^ + tp^y) (33) 
(34) 
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where P is causal and N is even, and both are real. Computing the derivatives we get 

^P n = P even + *N (35) 

j^Pu = Podd - iN (36) 

^P2l = -Podd-^~N (37) 

^P 22 = -P even + iN (38) 

where 

Peven (x, y) = ^ [P (x, y) + P (y, x)} ; P odd (x, y) = - [P (x, y) - P (y, x)} (39) 
Observe that P ret = NP (x,y). 

3 Translation-invariant solutions 

Translation invariant solutions are functions only of the relative variable x — y and may be 
Fourier transformed 

G ah (x - y) = [ e**"') G a6 (p) (40) 

The Fourier transform of an even (odd) kernel is an even (odd) function of p. If a kernel is 
real, the real (imaginary) part of its Fourier transform is even (odd). Vice versa, if a kernel is 
imaginary, then the real (imaginary) part of its transform is odd (even). 

It follows that, since Peven (x — y) is real and even, P eV en (p) is also real and even, while 
since P dd (x — y) is real and odd, P dd (p) is odd and imaginary. We may write 

Podd (p) = «7T7 (p) sign (41) 

therefore 7 (p) is real and even. 
3.1 The density of states 

Let us introduce the density of states A (p) out of the Fourier transform of the Jordan propa- 
gator G = G 21 - G 12 



G(p) = 27rA(p) sign p u (42) 
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The Jordan and retarded propagators are related through G (p) = 2 lmG ret (p) . From Eq. (|3T|) 



G 



ret 



:-i) 



W + 77 p - (p)1 



(43) 



Because of the retarded boundary conditions on G re t, it is understood that p° is replaced by 
p° + ze, e — > 0. We must distinguish two cases. As we shall see below, in the LO approximation, 
Zb = 1 and P ret (p) =0. In this case we have the explicit expression for G ret 



G ret (?) 



(p° + ie) + p 2 + x 



(LO) 



Therefore 



and 



G(p) =2 lmG ret (p) = 2tt sign (p°) 5 (p 2 + X ) (LO) 



A(p)=5(p 2 + x) (LO) 
In all higher approximations, P ret (p) ^ 0. We get 



G(p)=2\ G ret (p) | 2 ImP odd (p) (iVLO and higher) 



that is 



(44) 



(45) 



(46) 



(47) 



A (p) = \G ret (p) | 2 7 (p) (iVLO and higher) 



(48) 



3.2 The distribution function and the fluctuation-dissipation rela- 
tion 



Now consider the (12) component of Eq. (150) 



-^D (p) G 12 (p) + i [P n (p) G 12 (p) + P 12 (p) G 22 (p) 

Z_b iv 



Introducing the advanced propagator 
rewrite this as 



(49) 
i) (G 22 - G 12 ) (G adv (p) = [G ret (p)}*) we may 



[Gret (PT 1 = ^Pl2 (P) G adv (p) =► G 12 (p) = ^P 12 (p) |G ret (p)| 2 



and transform this into 



G 12 (p) = [-7T7 (p) sign (p°) + N (p) 



AJp) 

t(p) 



(50) 



(51) 
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In other words, for a translation invariant solution we must have 

G 12 (p) = 2tt F 12 (p) A (p) 
where (recall that sign (p°) = 1 — 26 {—p )) 

F 12 (p) = (-p°) + / (p) 
Comparing both expressions for G 12 we get 

~N(p) 



f(p) 



7T7 (p) 



- 1 



(52) 



(53) 



(54) 



It is more common to write this as 



N(p)=7r 7 (p)[l + 2/(p)] (55) 

whereby we recognize the fluctuation - dissipation relation. 

Given the Jordan and negative frequency propagators, it is easy to find all the others. 
In particular, the positive frequency propagator G 21 (p) = G + G 12 = 2ir F 21 (p) A (p) , with 
F 21 (p) = 9(p°) + f(p). 

3.3 A necessary condition for translation invariant solutions 

The expression for / (p) above (eq. (|54])) is equivalent to the identity 

P12 (p) G 21 (p) - P n (p) G 12 (p) = (56) 
I]) and (|54]) we get 



Indeed, from eqs. ( |3q 



A 2 (P) G 21 (p) - P 21 (p) G 12 (p) = P 12 (p) (9 (p°) - P 21 (p) (-p°) + [Pia (p) - (p)] / (p) (57) 
but also 



P12 (p) - P21 (p) = 2NP odd (p) = 2*A% 7 (p) sign (p°) 
P12 (p) (p°) - P21 (p) (-P°) = AT [P ocM (p) - *N (p) sign (p°) 

r n (p) 1 / 

= iNirj (p) 1 — sign (p° 

7r^/ In) V 



(58) 



7T7 (p) 

-2iA^7T7 (p) / (p) sign (p°) 



(59) 



Substituting these identities in Eq. (^) we get fl5*6"|), which is therefore a necessary condition 
for translation invariant propagators. 
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4 The large N approximation 



So far we have shown that translation invariant solutions are defined by the density of states 
A (p) and the distribution function / (p). To show that they correspond to thermal propagators, 
we must show that / (p) is necessarily of the form of a Bose - Einstein distribution 

f( p ) = Jel/fcl - lp (60) 

To do this, we need explicit expressions for the P kernels, which we can only find perturbatively. 

We shall adopt the large N approximation, which consists of taking the limit of iV — > oo 
in Yq and retaining only terms scaling like N (LO), 1 (NLO), N~ l (NNLO), etc. The key 
observation is that in any given Feynman graph each vertex contributes a power of N, each 
internal line a power of iV -1 , and each trace over group indices another power of N. We have 
both G and H internal lines, but the G lines only appear in closed loops. On each loop, the 
number of vertices equals the number of G lines, so there only remains one power of N from 
the single trace over group labels. Therefore the overall power of the graph is the number of 
G loops minus the number of H lines. Now, since we only consider 2PI graphs, there is a 
minimun number of H lines for a given number of G loops. For example, if there are two G 
loops, they must be connected by no less than 3 H lines, and so this graph cannot be higher 
than NNLO. A graph with 3 G loops can not have less than 5 H lines, and so on. 



4.1 The leading order approximation 



We conclude that Yq vanishes at LO, and therefore P a b 
the equations we need to solve become 



Q ao = 0. Under the ansatz Z B = 1, 



c ab D (p) - i 



G~ 



(p) = 



^{x-Ml}- h - f^- d G^(p)=0 



(61) 



(62) 



We disregard the auxiliary field propagator H ab , since to this order it is decoupled from the 
background auxiliary field and the other propagators. 

For the retarded propagator, we have the expression Eq. flUD, leading to Eq. Q35D for the 
Jordan propagator and ( 46|) for the density of states. Write Eq. (|6lD 



as 



D (p) G ab (p) = ic ab (63) 

Setting a = 2, b = 1, we see that G 21 can only be nonzero at the zeroes of D (p) = — p 2 — x, so we 
still can write G 21 (p) = 2vr F 21 (p) A (p) . Also G 12 (p) = G 21 (-p), and G 2 \p) -G 12 (p) = G (p) , 
so from Eq. ( fl5[) we conclude that F 21 (p) — F 21 (— p) = sign (p°) . Therefore we may write 
F 21 (p) = 9 (p°) + / (p) , with / (p) a real and even function. From these results, we may 
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write all propagators in terms of the distribution function / (p) . In particular, the Feynman 
propagator becomes 



G 



11 



-i) G ret + G 
(-1 



12 



(P 



,0 



is) +]? + x 
i-i) 



2tt 



(-P°) + f (p)l 5 (p 2 + x 



+ 2nf(p) 5(p 2 + x) 



_p02 _|_ ^2 _|_ ^ _ ^ £ 

So, assuming <i = 4 — e dimensions, we may evaluate 

^ d G l \p)=^ [M 2 +M 2 
where n is some (so far) arbitrary renormalization scale, 



(64) 



(65) 



2^r 



V 



2vre 



47T/i 2 



M 2 



d d p 

(2^r 



2tt/(p) <*(p 2 + x) 



(66) 
(67) 

A fl l~ -J 2- L 'J - (68) 

We are now confronted with the formal need to show that Eq. (|B"5D admits finite solutions 



The gap equation becomes 



My + Mj 







for x when e — > |54j] , as well as the physical need to show that the theory is reasonably stable 
against changes in the distribution functions / (p) [^5], |56], [57|, |58f . Let us interpret this equation 
as defining x as a function of Mj. Taking one derivative, we obtain 



A 



B 



X 



167r 2 e \Airn 2 



dx 



dM] 2^ 







(69) 



This suggests defining a background field and renormalization scale dependent effective coupling 
constant A from 



dx 
dMj 



(70) 



In other words 
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Now the gap equation 



-1 



Ml 



1 e 
+ 



A 2A 



B 



X 



shows that x — when Mj = M^ rit , given by 

l Ab 

So assuming M^ rit to be finite, we may rewrite the gap equation as 



1 e 
- + 



A 2A 



B 







(72) 



(73) 



(74) 



We see that it is possible to find a solution with finite propagators and x f° r an y distribution 
function /, provided Mj is finite. On the other hand, (weak) thermalization would require that 
any solution eventually converges to the thermal form Eq. (^Dj). Therefore, we conclude that 
the LO system does not thermalize. 

4.2 The NLO approximation 

Since the LO approximation admits a plurality of translation invariant solutions, it is necessary 
to go at least to NLO to study the issue of thermalization. There is only one NLO graph, 
consisting of a single G loop and a single H line (see Fig. 1). This graph leads to 



T^ LO = [ -i) (j^) 2 2N (V) Ca bc c def J d d xd d y H ad (x, y) G be (x, y) G c ? (x, y) 



h*o ( — 

Therefore, from 



P f \ 2 5T 

Pab {X, V) = T 



we get 



h5G ab (x,y)' 



Qab [x,y) = t- 



h5H ab (x,y) 



Pab (x, y) = -^c acd c bef H ce (x, y) G df (x, y) 



%h 

Qab {X, y) = ^TpCacdCbefG 06 (x, y) G df (x, y) 



(75) 
(76) 

(77) 
(78) 
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Since Q a b does not depend on H ab , we may solve the corresponding equation 



1 



A 



B 



—c ah 5(x,y) - i H 1 (x, y) + Q ab (x, y) = 

' ab 



in closed form. First, let us Fourier transform 

Ih- 1 



* {p) ~ { H g 2 i t^ + q 2 



(79) 



(80) 



To lowest order in 1/N, we may use the LO propagators to compute the Q's. In particular, we 
get 



where 



Qvu (p) 



Qn — A 4 e Qvu + Q/ii + Q/n 



-ift\ , /" o^o 



(2n) d (p - qf + x - <? 2 + X ~ ie 

d d q 2nf(q) 6 (q 2 + g 
(27r) d (p - g) 2 + x - 



=^i^f 4?r2 / ^) 5 2 +*) ftp-*) 6 o - + 

We assume that Q/n and Qfli are well defined, and compute 



X 



Qvn (p) 



ft r 



i + 



x ^ fx(l — x)p 2 + x\ 



167T 2 e JO \ 47T/i 2 

Recalling the renormalization condition Eq. (|7T|), we find 



(81) 

(82) 
(83) 
(84) 



(85) 



-!- + fi e Q V ii 
Xb 



i ftr 



- + 



i + 



x 



A lQn 2 e \4ir 



= -T + Qvrenll 



dx 



l+x{l-x) l^r 



-6/2 



(86) 



Qvrenii is explicitly finite. Let us call Q re „n = Qwenii + At e Q/n + Qfn • Since Q22 = —Qn, 
we may also write 



-1 , _ -1 , _ 

"7 r V22 — — : r Qren22 

Ab A 



i7) 
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where Q ren 22 = -Qv ren u + ^ e \~Qfii + Qfii\ ■ Introduce Q ret = Q ren ii + Q12, Qadv = 

Qrenll + QlX = Qret an d Write tnis aS (recall that Qren22 + Qrenll = Q22 + Qll = ~ (Ql2 + Q21)) 



H- 1 



ab 



(P) 



-I) 



T + Qret — Ql2 



Q 



12 



Q 



21 



(x + Qadv + Ql2 



Call 



(88) 



s- 1 (p) 



det 
1 



H- 1 

■)- Qret 



ab 



(P) 



A 



Qret — Ql2 



+ Qadv + Ql2 



+ Q12Q2I 



X 

"I - Qret 



+ Q 



V ~)~ Qadv 



+ Ql2 [Qret — Qadv] + Ql2 [Q2I _ Q 



12 



Observe that S (p) is real, finite and positive definite. Then 

' 1 + Qadv + Ql2 Ql2 



H ab (p) = z'E (p) 



Q 



21 



— + Qret — Q 



12 



?9) 



(90) 



4.3 The NLO density of states 

Beyond this point, our analysis will not depend upon the details of the NLO approximation, 
but only on a few structural features. One of these features is the fact that the NLO density 
of states is nonvanishing for — p 2 > 9x- To establish this fact, it is enough to look at P2i(p). 
Recall that, from Eq. (P7|), we know that P2i{p) is imaginary. Its odd part determines the 
kernel P dd, and its even part the noise kernel N. P dd determines 7 (p) through Eq. (|4T|). 
For an actual translation invariant solution, N and 7 are related by the fluctuation dissipation 
relation. The density of states A and 7 are related through Eq. (pE8p. It is clear from this 
equation that their zeroes are exactly the same, and so we only need to show that 7 is non 
vanishing. 
At NLO 



^ ^ ir/(^ 2IWG21(p - 9) 

= iS^k^ (q)an(T)an(q - r)G ' n(r - q) (91) 

To simplify the analysis, we can make the rather drastic approximation S (q) ~ A 2 , and use the 
LO propagators. 
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ih 2 X 2 r d d qA (q) d d rA (r) d d sA (s 



P21 ip) 



2Z% 



(27r) d_1 (27r) d-1 (2?r) 



— 5(g + r + s-p) 



* (?°) + / (9)] [O (r°) + / (r)] [0 + / (a) 



It is convenient to write the integrals in terms of future oriented momenta only. We get 

P21 (p) = z ^r~ [ J 3 + 3/ 2 + 3/i + Jo] 



2Z 



B 



where 



h = j DqDrDs 6 (q + r + s - p) [1 + / (q)} [1 + / (r)] [1 + / (s)} 
h = j DqDrDs S (q + r - s - p) [l + f(q)) [l + /(r)]/(s) 
/! = y DqDrDs8(s-q-r-p) [1 + / (s)] / (r) / (g) 



/ = | DqDrDs 6 {q + r + s +p) f {q) f (r) f {s) 



and we have defined 



(2vr) c 



T ^ (<?°) A (g) 



(92) 

(93) 

(94) 
(95) 
(96) 
(97) 

(98) 



It is clear that each of these integrals is nonnegative, so we only must show that for arbitrary 
p at least one is nonzero. 

Let us assume the LO density of states within the integrand, so that momenta q, r and s 
are on-shell. Assume — p 2 > 9%. If p° > 0, 1% and I2 are nonvanishing, while I\ and I2 are zero. 
If p° < 0, it is the other way round. If p° > 0, moreover, J 3 (p) > J (— p) and I 2 (p) > Ii (— p), 
and so both the odd and even parts of P\ 2 are nonvanishing, as we wanted to show. The fact 
that they are not only nonzero but actually proportional to each other only obtains for a special 
form of / (p) , indeed, a thermal form. We shall show this in next Section. 

It must be observed that in going from P\ 2 to P dd there is an extra factor oll/N involved 
(cfr. Eq. ([36])), and so the off-shell density of states, while non zero, is of higher order in 1/N. 



5 The only translation invariant solutions to NLO are 
thermal 

We may now show that the only translation invariant solutions of the NLO equations are 
thermal. The solutions must satisfy the identity Eq. (|56|), which becomes 
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/ 4\ d [^ 12 {q) ° 12 °° ~ q) ° 21 {p) " R21 {q) ° 21 °° ~ q) ° 12 = (99) 



(2*) 

^From the explicit solution for H 

d d q 



J — ^ £ (g) [Q 12 (g) G 12 (p - q) G 21 (p) - Q 21 (g) G 21 (p - g) G 12 (p)] = (100) 



(27T) 

Finally, use the iVLO approximation to Qafe 



/ 7^7^ S (?) I Gl2 W Gl2 - r ) Gl2 ( p " G21 ^ " G21 ^ G21 ( g - r ) ^) G21 ( p " g ) ° 12 ® 



(101) 



It is more usual to write this as 



= W / ( f q f *"„ f ' £ (p - g) 5 (g + r + 3 - p) {G 12 (g) G 12 (r) G 12 ( a ) G 21 (p) 
7 (^tt) l^ 71 ") l^ 71 ") 

-G 21 (g)G 21 (r)(g)G 21 (,)G 12 (p)} (102) 

We recognize the usual Boltzmann collision term, with £ (p — g) playing the role of cross section. 
Thus the only solutions must be thermal. 

There is one important observation to be made. Since this term is of order 1/N, in a strict 
power expansion we would use the LO density of states Eq. ( |46|) in the propagators. In the 
Introduction, we described this as making a "true" approximation . In practice, it means that 
we would put all momenta on mass - shell. Then only binary collisions would be possible, and 
the Boltzmann term would admit solutions with nonvanishing chemical potential. These do not 
exist in the exact theory, and so the "true" approximation does not describe thermalization. 

We see, however, that if we keep the NLO density of states the problem disappears. In 
the Introduction, we called this procedure the "exact" way. Because the density of states 
is nonzero everywhere, it is possible for one on - shell particle to decay into three off - shell 
ones, or vice versa, for one off-shell particle with momentum — p 2 > 9% to decay onto three 
on-shell particles (to show that this possibility is indeed open we have shown explicitly in the 
last Section that the NLO density of states is nonvanishing in this region). Particle number is 
no longer conserved, and only zero chemical potential is allowed. 

Thermalization is described, but only as a higher than NLO phenomenon, since it depends 
on NLO corrections to the density of states within an expression which is itself a NLO con- 
struct. We must consider if some of the NNLO terms we have left out may not bear on this 
process at a similar level. We shall see below that this is indeed the case. 
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6 Relaxation of the chemical potential 



We have seen so far that if nonperturbative corrections to the density of states are allowed, then 
the only translation invariant solutions to the Schwinger-Dyson equations are thermal propa- 
gators with vanishing chemical potentials. These are also the only true equilibrium solutions 
of the full theory. However, since we have stepped beyond the strict NLO approximation, we 
must ask ourselves if terms we have discarded would not affect thermalization at a similar level. 
To investigate this question, we shall consider the relaxation of the chemical potential. 

To place the issue in the simplest possible context, we shall assume that only long wave- 
lengths are involved, and adopt the Kadanoff-Baym approach. That is, we shall assume that 
all two-point functions G (x, y) may be written as functions of a relative variable u = x — y and 
a center of mass variable X = (x + y) /2, and that the X dependence is weak, loosely meaning 
that d x G ~ (l/N)d u G. 

We may introduce a mixed representation G (X, p) by performing a Fourier transform on 
the u variable (cfr. Eq. (|40|)). For fixed X, the manipulations in Section III are still valid, only 
now both the density of states and the distribution function display an extra X dependence 
(we refer the reader to ref. J2(| and references therein for a detailed discussion). 

To obtain the dynamics of the distribution function / (X, p) , write the Schwinger-Dyson 
Eq. (H) for G 21 (x, y) as 



±-D (x, y) G 21 (x, V)-^J d d z [P 21 (x, z) G n (z, y) + P 22 (x, z) G 21 (z, y)] = (103) 

Recall that D(x,y) = [d 2 — x(x)]5(x,y) and Fourier transform with respect to the relative 
variable in each case to get, up to 1/N terms 



P 21 (X,p) G (X,p) + P 22 (X,p) G (X,p) 



1-n + iL] G 21 [X,p) = -1 

where Q = p 2 + x {X) , and L is the Vlasov operator 

d 1 dx d 
~ P dX ~ 2dXdp 
Separating the imaginary part, we get the Boltzmann equation 



(104) 



(105) 



Lf = l 



col 



where (recall eqs. (j35| ) to (p^)) 

-iZ B 
2N 



l col 



P 21 (X,p) G (X,p) - P 12 (X,p) G (X,p) 



(106) 



(107) 



Observe that different approximations yield different collision integrals. Of course, the 
condition I co i = for a translation invariant solution is just the necessary condition Eq. (|56|). 
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Manipulating the collision term as in the last Section, we get, to NLO 



(2Ti) d \ 2 h 2 f d d q d d r d d s , 

{G 12 (q) G 12 (r) G 12 (s) G 21 (p) - G 21 (q) G 21 (r) (q) G 21 (s) G 12 (p)} (108) 



where we leave the X dependence implicit, and have approximated £ ~ A 2 for simplicity. 

Since we are only interested in the relaxation of the chemical potential, we may linearize 
the Boltzmann equation around equilibrium. Write the distribution function as 

/ (X, p) = /„ (p) [1 + (1 + /o) 5f (X, p)] (109) 

where /o is a thermal distribution with vanishing chemical potential. Neglecting the variation 
in the density of states, which gives a higher order contribution, we get 

G 21 ' 12 = 2vrA (p) [9 (±p°) + f + f (1 + f ) Sf (X,p)] (110) 

A nonzero chemical potential corresponds to a p-independent perturbation 5f(X,p) = 
8fi (X). To isolate the dynamics of the chemical potential, we may integrate both sides of the 
Vlasov-Boltzmann equation with respect to p. Since the distribution function is even, there is 
no loss in resticting the integration region to p° > 0. We also write the collision integral in 
terms of integrals over positive energy momenta, to get 

fl s*w-^S' (m, 



where (recall the momentum space measure Eq. fl98|)) 

B = j Dpf (l + f )p° (112) 

J = J DpDqDrDs 5 (q — r — s — p) 

{/ (<?) [1 + / (r)] [1 + / (s)] [1 + / (p)) - [1 + / (q)] f (r) / (s) f (t)} (113) 

Linearizing the J integral we get 

J = -2K5fi(X) (114) 
K = J DpDqDrDs 5 (q - r - s - p) [1 + / (q)} [1 + / (r)] [1 + / (s)] [1 + / (p)] e^ 9 (115) 
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The integral is nonzero, provided we use the NLO density of states. If we used the LO density 
of states, then the argument of the delta function would never vanish. This means we predict 
chemical potential relaxation on a time scale which grows at least as N 2 , since to the 1/N 
factor in Eq. (|111|) we must add at least one more 1/N factor coming from the off-shell density 
of states. 

The problem is that there are terms, coming from the NNLO approximation to the 2PI 
CTP EA, that modify the collision term in such a way to contribute to chemical potential 
relaxation at the same accuracy. Therefore the NLO prediction for the relaxation time is not 
accurate. We shall conclude this Section by showing this explicitly. 



6.1 The NNLO approximation and particle number violation 

As we already remarked, the power of N associated to a given vacuum bubble is the number of 
G loops minus the number of H lines. Therefore, in the NNLO approximation we must look 
for graphs with one more H line than G loops. However, a 2PI graph with 3 G loops must 
have no less than 5 H lines, so the allowed number of G loops is only 1 or 2. There are only 
two graphs with 1 G loop and 2 H lines, but only one of them is 2PI. (see Fig. 2) Similarly, 
there is only one IP I graph with 2 G loops and 3 H lines (see Fig. 3). We therefore have two 
new graphs contributing to the IP I CTP EA. 

Variation of these graphs with respect to G yields two new contributions to the P a b (Figs. 
4 and 5) 

Let us consider Fig. 4, and replace the H lines by their expansion in powers of A (Fig. 6) 
The lowest order contribution yields the setting sun diagram (Fig. 7). 

Introducing the first correction to one of the H lines gives the graph in Fig. 8. Correcting 
both H lines gives the graph in Fig. 9. 

The second new graph in the 2PI CTP EA (Fig. 3) yields, upon variation, the graph in 
Fig. 5. 

Replacing the H lines by their lowest order expression, we find another contribution of the 
form of Fig. 8. The next order yields contributions proportional to Figs. 9 and 10. 

It has been shown in ref. (|2(J) that the graphs in Fig. 9 and 10, when translated in terms of 



the collision integral, describe scattering of 2 into 4 particles and vice versa. These scattering 
processes do not conserve particle number and therefore contribute to the relaxation of the 
chemical potential (we have shown explicitly the linearized collision operator in ref. (@))- 
The resulting contribution is at least of the same order of magnitude as that found by allowing 
a nonperturbative density of states in the NLO collision term. 

Therefore we conclude that the NLO prediction for the relaxation time of the chemical 
potential is not accurate. It is nevertheless remarkable that the NLO succeeds in predicting 



relaxation, in agreement with the claims of refs. [21| and 
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7 How does macroscopic irreversibility appear? 



As we mentioned at the beginning of this paper, the issue of thermalization in relativistic 
quantum fields, besides its concrete application to high energy plasmas, is relevant to the 
larger question of the origin of the thermodynamic arrow of time in physics. The LO Large N 
approximation does not break time reversal invariance (it does not lead to thermalization either) 
and one could jump to the conclusion that this is a general feature of the Large N perturbative 
scheme. However, as we have seen, the NLO or, at worst, the NNLO approximations show 
thermalization. It behooves us to identify at which point time-reversal invariance has been 
broken, and how. 

To identify the crucial assumptions, let us return to the form Eq. ( |T0D of the action for the 
theory. A variation with respect to the field yields the Heisenberg equations of motion for the 
field operator 



Multiplying from the left (say) by tp a (y), taking expectation values and summing over a, we 
obtain (recall eqs. flT6|) and ([24])) 



This Schwinger-Dyson equation is, of course, not a closed equation for the propagator; 
it only relates the propagator to a higher correlation function. If we wish to say something 
about this new correlation, one possibility is to repeat the argument. We may begin from the 
Heisenberg equation for the auxiliary field 



<9V (z) - X (x) ¥ a {?) = 



(116) 



hdlG 21 {x,y)-{ X {x)y a (x) <p a (y)) = 



(117) 



X (x) = Ml + -fy° (x) <p a (x) 



(118) 



and multiply by ip a (x) (p a (y) to get 



(X (x) V a (x) y\y)) = MlUG 21 (x, y) + ^L ^ {x) ^ (x)(p a {x) 



(119) 



or else we go back to Eq. (|116|) to get 




(120) 



Comparing Eq. ( |117|) to Eq. (p6|) , we see that 



1 



( X (x)p a (x) V a (y))=x(x)G 21 (x,y) 



1 



/ 



d 4 z [P\ (x, z) G u {z, y) + P\ {x, z) G 21 (z, y) 

(121) 



h 



N 



Using Eq. (p8|) for Xi an d comparing to Eq. (|119|) , we get 
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^| (x) <p\x) V a (x) ip a {y) 



^G 21 (x,x)G 21 (x,y) 



~N 1 ** 



P\ (x, z) G (z, y) + P\ (x, z) G (z, y^.22) 



To be more concrete, observe that if Wick's theorem held, then 

(<p b (x) y b {x)tp a (x) <p a (y)) ~h 2 (l + A) G 21 (x, x)G 21 (x, y) (123) 
This suggests writing 

b (x) ip\x)<f (x) V a (y)) = h 2 [G 21 (x, x)G 21 (x, y) + 1(7 (x, y)] (124) 

The NLO approximation consists in closing the Schwinger- Dyson hierarchy by writing the 
Pat's as in Eq. (|77|), whereby 



^-C (x, y) = -ij d 4 z [H 21 (x, z) G 21 (x, z) G 11 (z, y) - H 22 (x, z) G 22 (x, z) G 21 (z, y)] (125) 

Observe that the integrand vanishes when z° >z x°, y°. At this point, time reversal symmetry 
has been broken. 

Let us investigate further the mechanism for breaking time symmetry. First decompose the 
H propagators in singular and regular parts 



H ab (x, y) = H (x) c ab 5 (x - y) + Hf eg (x, y) 



Then C (x, y) is split into a reducible and an irreducible term 



C(x,y) 



-2i 



A 



H (x) G (x, x)G 21 (x, y) + C irr (x, y) 



(126) 



(127) 



B 



Time symmetry is broken because the irreducible term Ci rr (x, y) vanishes in the distant past 
while remains non-zero in the far future. 

For completeness, let us observe that Eq. fl2? ) shows that H (x) = i\b + O (X B ) , and 
therefore Eq. ( |127| ) leads to Eq. ( |123[ ) to leading order. 



In summary, the scheme works because at a crucial point some higher correlation (C in 
eq. ( |124|) ) is replaced by a perturbative expansion in terms of propagators. The replacement 
assumes that the irreducible part of the higher correlation (Cj rr in eq. ( |127|) ) vanishes in the 
distant past, in effect enforcing a variant of Boltzmann's molecular chaos condition. 

In order to restore time reversal symmetry, we ought to treat C (x, y) as a dynamical vari- 
able in its own right, for example, by defining a higher generating functional with a new non 



local source coupled to four fields. We refer the reader to ref. pO] for a fuller discussion of this 
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issue. 



Acknowledgements We acknowledge discussions of NLO scheme with Emil Mottola and 
Paul Anderson. EC is supported in part by CONICET, UBA, Fundacion Antorchas and AN- 
PCyT. BLH is supported in part by NSF grant PHY98-00967 and their collaboration is sup- 
ported in part by NSF grant INT95-09847. 

References 

[1] E. Shuryak, Two-stage equilibration in high energy heavy ion collisions, Phys. Rev. Lett 
68, 3270 (1992) 

[2] E. Shuryak and L. Xiong, Dilepton and photon production in the "hot glue" scenario, 
Phys. Rev. Lett 70, 2241 (1993). 

[3] T. Hirano, Is early thermalization achieved only near midrapidity in Au-Au collisions at 
v /s^7 = 130 GeV?, Phys. Rev. C65, 011901(R) (2001). 

[4] U. Heinz and P. Kolb, Early thermalization at RHIC, Nucl. Phys. A702, 269 (2002). 

[5] U. Heinz and S. Wong, Elliptic flow from a transversally thermalized fireball, |hep- 
| ph/0205055 . 

[6] R. Baier, A. Mueller, D. Schiff and D. Son, Bottom-up thermalization in heavy ion colli- 
sions, Phys. Lett. B502, 51 (2001). 

[7] R. Baier, A. Mueller, D. Schiff and D. Son, Does parton saturation at high density explain 
hadron multiplicities at RHIC?, [hep-ph/02042TT| . 

[8] S. Coleman and R. Jackiw and H. D. Politzer, Spontaneous symmetry breaking in the 
O(N) model for large N, Phys. Rev. D10, 2491 (1974) 

[9] R. Root, Effective potential for O(N) model to order 1/N, Phys. Rev. D10, 3322 (1974) 

[10] F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz and P. Anderson, Nonequilibrium 
quantum fields in the large-N expansion, Physical Review 50, 2848 (1994) 

[11] F. Cooper, Y. Kluger, E. Mottola and J. P. Paz, Quantum evolution of disoriented chiral 
condensates, Physical Review 51, 2377 (1995) 

[12] J. Cornwall, R. Jackiw and E. Tomboulis, Effective action for composite operators, Phys. 
Rev. D10, 2428 (1974) 

[13] J. Schwinger, Brownian Motion of a Quantum Oscillator, J. Math. Phys. 2, 407 (1961). 



25 



[14] L. V. Keldysh, Diagram technique for nonequilbrium Processes, Zh. Eksp. Teor. Fiz. 47, 
1515 (1964) [Sov. Phys. JETP 20, 1018 (1965)] 



[15 
[16 

[17; 

[18 
[19 
[20 
[21 
[22 
[23 



P. M. Bakshi and K. T. Mahanthappa, J. Math. Phys. (N.Y.) 4, 1 (1963), ibid, 4, 12 (1963) 

K. Chou and Z. Su and B. Hao and L. Yu, Equilibrium and nonequilibrium formalisms 
made unified., Phys. Rep. 118, 1 (1985). 

E. Calzetta and B. L. Hu, Nonequilibrium quantum fields, Phys. Rev. D37, 2878 (1988) 

L. Kadanoff y G. Baym, Quantum Statistical Mechanics (Addison - Wesley, New York, 
1962). 

P. Danielewicz, Quantum Theory of Nonequilibrium Processes I, Ann. Phys. (NY) 152, 
239 (1984). 

E. Calzetta, B-L. Hu and S. Ramsey, Hydrodynamic transport functions from quantum 
kinetic field theory, Phys.Rev. D61, 125013 (2000) 

J. Berges and J. Cox, Thermalization of Quantum Fields from Time-Reversal Invariant 
Evolution Equations, Phys. Lett. B517, 369 (2001). 

J. Berges, Controlled nonperturbative dynamics of quantum fields out of equilibrium, Nucl. 
Phys. A699, 847 (2002). 



L. Cugliandolo, D. Grempel, G. Lozano, H. Lozza and C. da Silva Santos, |cond- 



mat/01 10497 . 



[24] S. Jeon and L. Yaffe, From quantum field theory to hydrodynamics: transport coefficients 
and effective kinetic theory, Physical Review D53, 5799 (1996) 

[25] P. Arnold and L. Yaffe, Effective theories for real-time correlations in hot plasmas, Physical 
Review D57, 1178 (1998) 

[26] P. Arnold, D. Son and L. Yaffe, Effective dynamics of hot, soft non-abelian gauge fields: 
Color conductivity and log (1/alpha) effects, Physical Review D59, 105020 (1999) 

[27] P. Arnold, D. Son and L. Yaffe, Longitudinal subtleties in diffusive Langevin equations for 
non-abelian plasmas, Physical Review D60, 025007 (1999). 

[28] D. Bodeker, From hard thermal loops to Langevin dynamics, Nuclear Physics B559, 502 
(1999) 

[29] J. P. Blaizot and E. Iancu, The quark-gluon plasma: collective dynamics and hard thermal 
loops, Phys. Rep. 359, 355 (2002) 



26 



[30] E. Calzetta and B. L. Hu, Stochastic dynamics of correlations in quantum field theory, 
Phys.Rev. D61, 025012 (2000) 

[31] D. Litim and C. Manuel, Semi-classical transport theory for non-Abelian plasmas, Phys. 
Rept. 364, 451 (2002). 

[32] S. Mrowczynski, Transport theory of massless fields, Phys. Rev. D56, 2265 (1997). 

[33] J. Hartle and G. Horowitz, Ground-state expectation value of the metic in the 1/N or 
semiclassical approximation to quantum gravity, Phys. Rev. D24, 257 (1981). 

[34] O. Eboli, R. Jackiw and S-Y. Pi, Quantum fields out of thermal equilibrium, Phys. Rev. 
D37, 3557 (1988) 

[35] F. D. Mazzitelli and J. P. Paz, Gaussian and 1/N approximation in semiclassical cosmology, 
Phys. Rev. D39, 2234 (1989) 

[36] S. Habib and Y. Kluger and E. Mottola and J. Paz, Dissipation and Decoherence in Mean 
Field Theory, Phys. Rev. Lett. 76, 4660 (1996) 

[37] D. Boyanovsky, H. de Vega, R. Holman, S. Kumar and R. Pisarski, Non-equilibrium evo- 
lution of a 'tsunami': dynamical symmetry breaking, Physical Review D57, 3653 (1998) 

[38] D. Boyanovsky, H. J. de Vega, R. Holman and J. F. J. Salgado, Nonequilibrium Bose- 
Einstein condensates, dynamical scaling, and symmetric evolution in the large N $ 4 theory, 
Phys. Rev. D59, 125009 (1999). 

[39] D. Boyanovsky, H. J. de Vega, R. Holman and J. F. J. Salgado, Analytic and numerical 
study of Preheating dynamics, Phys. Rev. D54, 7570 (1996) 



[40] D. Boyanovsky, F. Cao and H. de Vega, Inflation from Tsunami waves, |astro-ph/0102474| . 



[41] S. A. Ramsey and B. L. Hu, O(N) Field theory in curved spacetime, Phys. Rev. D56, 661 
(1997) 

[42] S. A. Ramsey and B. L. Hu, Nonequilibrium inflaton dynamics and reheating, Phys. Rev. 
D56, 678 (1997) 

[43] P. Greene, L. Kovman, A. Linde and A. Starobinsky, Structure of resonance in preheating 
after inflation, Phys. Rev. D56, 6175 (1997). 

[44] L. Bettencourt and C. Wetterich, Time evolution of correlation functions in Non- 
equilibrium Field Theories, Phys. Lett B430, 140 (1998). 

[45] D. T. Son, Reheating and thermalization in a simple scalar model, Phys. Rev. D54, 3745 
(1996) 



27 



G. Aarts, D. Ahrensmeier, R. Baier, J. Berges and J. Serreau, Far-from-equilibrium dy- 
namics with broken symmetries from the 2PI — 1/N expansion, |hep-ph/0201308| . 

B. Mihaila, T. Athan, F. Cooper, J. Dawson and S. Habib, Exact and approximate dy- 
namics of the quantum mechanical O(N) model, Phys. Rev. D62, 125015 (2000). 

B. Mihaila, J. Dawson and F. Cooper, Resumming the large-A r approximation for time 
evolving quantum systems, Phys. Rev. D63, 96003 (2001). 

L. Bettencourt and C. Wetterich, Time evolution of correlation functions for classical and 
quantum anharmonic oscilators, hep-ph/9805360 v2. 



G. Aarts, G. Bonini and C. Wetterich, Exact and truncated dynamics in nonequilibrium 
field theory, Phys. Rev. D63, 25012 (2000). 

G. Aarts, G. Bonini and C. Wetterich, On thermalization in classical field theory, Nucl. 
Phys. B587, 403 (2000). 

K. Blagoev, F. Cooper, J. Dawson and B. Mihaila, Schwinger-Dyson approach to nonequi- 
librium classical field theory, Phys. Rev. D64, 125003 (2001). 

G. Aarts and J. Berges, Classical aspects of quantum fields far from equilibrium, Phys. 
Rev. Lett. 88, 41603 (2002) 

H. van Hees and J. Knoll, Renormalization in self-consistent approximation schemes at 
finite temperature: Theory, Phys. Rev. D65, 25010 (2001). 

M. van Eijck, D. O 'Connor and C. Stephens, Heating field theory the environmantally 
friendly way, [hep-th/9310199| vl 

D. O 'Connor and C. Stephens, Environmentally friendly renormalization, Int. J. Mod. 
Phys. A9, 2805 (1994) 



D. O 'Connor, C. Stephens and A. Bray, Dimensional crossover in the large N limit, cond- 



mat/9601146> l 



D. O 'Connor, C. Stephens and J. Santiago, Dimensional Crossover in the Non-linear sigma 
model, |cond-mat/0a02205|vl. 



28 



Figure 1: Only NLO contribution to the 2PI effective action. The full line denotes a G 
propagator, while the dotted line stands for H 




Figure 2: Only NNLO contribution to the 2PI effective action with a single G loop. 




Figure 3: Only NNLO contribution to the 2PI effective action with 2 G loops. 



29 



\ 

\ 



/ 

/ 



Figure 4: Contribution to P a b from the variation of Fig. 2. 



Figure 5: Contribution to P a b from the variation of Fig. 3. 




Figure 6: Expansion of H in powers of the coupling constant. 



Figure 7: Setting sun graph 
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Figure 8: First correction to Fig. 4 



Figure 9: Second correction to Fig. 4 



Figure 10: Higher order correction to Fig. 5 
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